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A hybrid model of the Deutsch-Jozsa algorithm is presented, inspired by the proposals of hybrid 
computation by S. Lloyd and P. van Loock et. al. The model is based on two observations made 
about both the discrete and continuous algorithms already available. First, the Fourier transform is 
a single-step operation in a continuous-variable (CV) setting. Additionally, any implementation of 
the oracle is nontrivial in both schemes. The steps of the computation are very similar to those in 
the CV algorithm, with the main difference being the way in which the qunats, or quantum units of 
analogic information, and the qubits interact in the oracle. Using both discrete and continuous states 
of light, linear devices, and photo-detection, an optical implementation of the oracle is proposed. 
For simplicity, infinitely squeezed states are used in the continuous register, whereas the optical 
qubit is encoded in the dual-rail logic of the KLM protocol. The initial assumption of ideal states 
£18 qunats will be dropped to study the effects of finite squeezing in the quality of the computation. 
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I. INTRODUCTION 

The model proposed by Knill, Laflamme and Milburn 
[U E] presents an efficient way to implement a quantum 
computer using optical qubits, single-photon sources, lin- 
ear devices and photo-detection, knovi^n as the KLM pro- 
tocol. This well established model uses nonlinearities 
hidden in the measurement process to produce the inter- 
actions needed between the optical modes, together with 
offline state preparation and teleportation to enhance the 
probability of success of the conditional gates of the pro- 
tocol, and error correction codes, achieving a scheme that 
is both fault-tolerant and robust against errors [31 13] . 

On the other hand, the unconditional nature in which 
entanglement with continuous variables can be obtained, 
and thus the unconditionalness of the operations per- 
formed over CV together with the need of just single- 
mode nonlinearities for universal quantum computation, 
are two of the most prominent reasons to adopt this set- 
ting over a discrete- variable one [51 Besides, the idea 
of studying quantum algorithms, generally designed for 
qubits, in a CV setting is at the same time challenging 
and of great importance for the development of quantum 
information theory. This generalization may bring into 
light new algorithms with a more natural formulation 
using CVs. For example, any algorithm based on the 
quantum Fourier transform, such as the Deutsch-Jozsa 
algorithm or Shor's algorithm, is easily performed in a 
continuous setting, where such operation is a single-step 
operation (7|. 

Focusing on the Deutsch-Jozsa algorithm, it can be 
seen that its original discrete formulation [8 requires the 
production of superpositions of several qubit states via 
a quantum Fourier transform, a task that needs a lot 



of resources and computational power, even for a small 
number of qubits [Hj- The translation of the algorithm 
into a CV setting seems more feasible and practical, due 
to the ability to produce the superpositions needed with 
a very low number of resources, nonetheless, how to im- 
plement the oracle in such continuous scheme is still an 
unresolved problem 1 10] . 

However, we can take advantage of both schemes by 
means of some sort of hybrid gates, i.e. a set of trans- 
formations that acts on both qubits and CV states [11) . 
In this paper we use a hybrid model of computation, in 
an all-optical setting, giving an explicit optical circuit 
for the action of the oracle on both registers involved in 
the computation. Such model is based on proposals by 
S. Lloyd [TT] and P. van Loock et. al. [H], but differs 
from them by utilizing only linear optical devices, single- 
photon sources, offline squeezing and photo-detectors. 
In the first part of the work, we use idealized infinitely 
squeezed states as the computational CV basis for sim- 
plicity in our calculations. By the end, we shall drop 
such restriction to discuss the effects of realistic Gaus- 
sian states (that is states with finite squeezing) through- 
out the computational process. 

The structure of the paper is as follows. First, in Sec- 
tion [n] we give a quick review of the original algorithm 
and its implementation in the continuous-variable case. 
Section |III| presents the algorithm in a hybrid and all- 
optical setting, paying special attention to the explicit 
implementation of the oracle. Section IV considers im- 



perfect CV states in the hybrid model of the algorithm, 
in order to analyze the possible distortions at the end 
of the computation. Finally, Section jV] closes the paper 
with some concluding remarks. 
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II. FUNDAMENTALS OF THE 
DEUTSCH-JOZSA ALGORITHM 

In this section we present a review of the basic aspects 
of the Deutsch-Jozsa problem according to its original 
formulation. Then, the explicit algorithm is given for 
CVs as presented in Ref. [lO] . 

A. Overview of the algorithm 

In the Deutsch-Jozsa problem [H] , we are given a black 
box quantum computer, known as the oracle, that imple- 
ments the function / : {0, 1}" — ?► {0, 1}. We are promised 
that the function is either constant (0 on all inputs or 1 
on all inputs) or balanced (returns 1 for half of the input 
domain and for the other half;) the task then is to de- 
termine if / is constant or balanced by using the oracle. 
Although it is of little practical use, it is one of the first 
examples of a quantum algorithm that is exponentially 
faster than any possible deterministic classical algorithm. 
It also provided inspiration for Shor's and Grover's algo- 
rithm [13l [H] , two of the most revolutionary quantum 
algorithms. 

For a conventional deterministic algorithm where n is 
the number of bits, 2"~^ -I- 1 evaluations of / will be 
required in the worst case. That is to say that in or- 
der to prove that / is constant, just over half the set 
of outputs must be evaluated and found to be identical 
(remembering that the function is guaranteed to be ei- 
ther balanced or constant, not somewhere in between.) 
The best case occurs when the function is balanced and 
the first two output values chosen happen to be different. 
For a conventional randomized algorithm, a constant k 
evaluations of the function suffices to produce the correct 
answer with a high probability (failing with probability 
£ < 1/2*''"^.) However, k = 2"~^ -I- 1 evaluations are 
still required if we want an answer that is always correct. 
The Deutsch-Jozsa algorithm produces an answer that is 
always correct with a single evaluation of /. 

B. The algorithm with CVs 

Implementing a superposition for a large number of 
qubits in a discrete setting, although possible, seems non- 
practical. This represents an important drawback for the 
implementation of many quantum algorithms that highly 
depend on such operation. The Deutsch-Jozsa algorithm 
is one of them and in the present section we show an 
alternative to apply the algorithm using qunats. In this 
case, the advantage of using a continuous setting comes 
from the low computational resources needed to imple- 
ment a continuous superposition; in fact, this is easily 
achieved through a Fourier gate on the incident qunat 
state [7]. 

In the CV case of the Deutsch-Jozsa algorithm [TU] . 
Alice and Bob play the following game: Alice chooses 



a random real number q between — oo and -I-cxd (or be- 
tween —L and +L in practice,) and then, she asks Bob to 
evaluate the function f{q) with possible outcomes and 
1. Notice that this game implies that, classically, Alice 
could only discover the nature of the function with total 
certainty if she asks Bob an infinite number of times! A 
quantum algorithm, however, only requires one evalua- 
tion of the function to determine if it is either constant 
or balanced, thus providing a significant speedup over 
any classical algorithm. 
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Figure 1. Quantum circuit for the Deutsch-Jozsa algorithm 
with continuous variables. 

To store her query, Alice has a qunat register to store 
her query in, and so does Bob, who stores the answer of 
the oracle in it. The quantum circuit representing the ac- 
tion of the algorithm on both qunat registers is shown in 
figure [T] Initially, the states are prepared in the ideal CV 
states |qo) and |7r/2) for Alice and Bob, respectively, both 
in the g-quadrature basis of the electromagnetic field. 
Afterwards, they both produce a superposition in each 
register using Fourier gates, producing the state 

^ j dgdg'e2**«+"«'|g)|(z')- (1) 

The oracle then performs the nondemolition operation 
between the qunat states \qi)\q2) ki) 192 + /(gi)), 
where the information about / is stored in the second 
variable. 

^ y d(? d(7'e2'®«+"«' \q) \q' + /(g)) 

= -^1 dge2^*«(-l)^'(^)|g)-Fk/2). (2) 

At this point, Bob's register has done its work and re- 
mains unchanged until the end of the computation. We 
apply one last (inverse) Fourier gate to Alice's register 
before the measurement stage 

IQ) = ^ / d(?dg'e2*'?(*-9')(~l)/(?) \q') . (3) 

Finally, Alice makes a projective measurement in her 
qunat state. She projects it onto the initial state |go), 
using the operator 

/•go+Ago/2 

nA,o = / dfc |fc) (fc| , (4) 
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which takes into account certain spread Ago in the mea- 
surement process. The state \Q) becomes then 

1 rqo+Aqo/2 /■ 

riA^o \Q)^- dfc / dge2^«(*-^-)(-l)^(«) |fc) . 

^ Jqo-Aq„/2 J 

(5) 

Let us first suppose the function is constant. In that 
case, /(g) = ±1 necessarily, so the projection in ^ can 
be simpHfied to 

nAgo IQ) = ± |go) , (6) 

and then Ahce measures with total certainty her initial 
quadrature state |go)- In the other case, when the func- 
tion is balanced, the output quadrature would have any 
value but qq, yielding a zero measurement. 

In practice, the preparation of nonideal computational 
states and finite precision in measurement pose severe 
problems for the experimental realization of this scheme, 
specially because the constant or balanced nature of the 
function would not be as straightforward to differentiate 
as in the idealized case presented above, the outcome of 
the measurement would not be as certain. Additionally, 
this model does not address the question of the imple- 
mentation of the oracle, but suggests that using some 
kind of hybrid setting (i.e. one that combines discrete 
and continuous states,) one may improve the algorithm. 

III. DEUTSCH-JOZSA ALGORITHM IN A 
HYBRID SETTING 

In the hybrid computer proposed in Ref. nec- 
essary and sufficient conditions are given to build up a 
universal set of gates operating in both a qubit and a CV 
state. Actually, the construction of such set is completed 
with just one interaction gate between both states, be- 
cause universality is achieved by the repeated application 
of such interaction gate and the single-mode gates of each 
setting. 
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Figure 2. Quantum circuit for the Deutsch-Jozsa algorithm 
in a hybrid setting. 

Let us consider the Deutsch-Jozsa problem in a hy- 
brid setting. We use a qunat state to store Alice's query 
and a qubit to store the result of Bob's evaluation. The 
quantum circuit is given in figure [2j We keep the contin- 
uous part to implement the Fourier gate, which is much 
more efficient than the discrete superposition achieved 



with consecutive Hadamard gates but we also use a 
discrete register in which it is easier to store the answer 
given by the oracle. 

A. Optical implementation of the algorithm 

The all-optical setting we present here starts with an 
ideal [T5] qunat state |go) in Alice's register, and a dual- 
rail qubit in Bob's register. The latter requires two paths 
or rails through which a single photon travels, a \0)q 
qubit is encoded if the lower path is populated |0, 1), and 
a qubit is obtained when the upper rail carries the 
photon [3] . For the particular case of this algorithm. Bob 
produces a qubit before the computation. 




Figure 3. Optical circuit of the hybrid Deutsch-Jozsa algo- 
rithm. The explicit implementation of the oracle is shown. 

Having prepared both registers, Alice and Bob pro- 
duce superpositions on each of them. For the CV state, 
she uses a phase shifter with 93 = 7r/2 to implement the 
Fourier gate Oil]. On the other hand, Bob performs the 
Hadamard gate using a 50-50 beam splitter {9 — n/A.) 
The optical circuit is shown explicitly in figure [3] The 
resulting state just before the oracle operation is then 

(-^ J dge2^'«^|g))®-^(|0),-|l),). (7) 

In this setting, we require that the action of the oracle 
transforms the state (TtI into 




in analogy to the action of the oracle in the CV setting. 
For now, it is all we are saying about the oracle, the 
complete analysis is given in the next section. 

As usual, Bob's register is left unchanged after the or- 
acle operation, and Alice performs one last superposition 
using an inverse Fourier gate (via a phase shifter with 
ip — — 7r/2.) Actually, these final steps are identical to 
those in the CV implementation, i.e. steps ^ through 
Once more, if the function is constant, Alice mea- 
sures her initial state |go) with total certainty. In con- 
trast, she measures zero if the function is balanced. 
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B. Operation of the Oracle 

Now, the crucial step of the algorithm occurs in the 
oracle, where the registers interact with each other. The 
hybrid nondemolition operation required in this point is 
given by 

\q)\x),^\q)\x®fiq))^, (9) 

where \q) is an infinitely squeezed state, and \x)q is an 
optical qubit. The operation x © /(g) means addition 
modulo 2. Now we make an important observation. As 
f{q) can only be or 1, this operation looks very much 
like a CNOT gate in the discrete setting. Let us use this 
intuition to construct the oracle with linear optical de- 
vices, photon sources and photo-detection, as shown in 
the boxed section of figure |3] 

After the superposition stage, the oracle evaluates the 
function once in the superposition state created by Alice. 
The result (0 or 1) is used to create another dual-rail 
control qubit, which is conditioned by the feed- forward 
given by the oracle: 

1. If f{q) — 0, we create a state |0)^ — |0, 1), that is, 
the photon source of the top rail does not produce 
anything, and the source of the bottom path creates 
a single photon, 

2. or if f{q} = 1, we produce a single photon in the 
upper arm and no photons in the other, obtaining 
the qubit = |1,0). 

This method allows the qunat to pass the oracle un- 
changed, the only operation made on it transmitted clas- 
sically to the ancillary control qubit. Therefore, the op- 
eration ([9]) becomes simply a CNOT gate between qubits 

\x)c\y)g^\^)c\y®^),, (10) 

with the information about / hidden in the value of x. 
The optical implementation of a CNOT gate can be effi- 
ciently performed using the KLM protocol of computa- 
tion together with error correction codes [U [5] . 

Remember that Bob's register is in the state (|0)^ — 

|l)^)/-\/2 before the oracle action. The CNOT gives us 
then the desired result: if f{q) — 0, Bob's qubit is left 
unchanged, whereas if f{q) = 1, it gains a sign. The net 
action is then a phase (— as expected. 

IV. EFFECTS OF GAUSSIAN STATES IN THE 
COMPUTATION 

Instead of using the unphysical, infinitely squeezed 
states, we replace them for the more realistic Gaussian 
states, that represent Gaussian fields such as the vacuum 
state |0) or states with finite squeezing iS(s)|0). The lat- 
ter are produced with nonlinear optical processes such as 
optical parametric amplification, which allows us to gen- 
erate squeezed vacuum states |16j . The wave function of 



a general Gaussian state is given by the distribution of 
variance s^/2 

Gs(g) = (gk)-(^s)-l/2e-''/^^ (11) 

and thus, the vacuum squeezed state \s) — S{s)\0) can 
be written as 

\s} = J dqGM\q)- (12) 

The parameter s is related to the squeezing parameter 
through r = ln(s), and thus we shall call it also the 
squeezing parameter. 

Now we can follow the steps of the algorithm with these 
realistic states in Alice's CV register. The initial state is 
then |s) (g) \l)q or 

(^(nsr'/'Jdqe-'^'/^'lq)'^®]!)^. (13) 

Both the continuous and discrete superpositions are per- 
formed in the respective register, transforming the state 
into 

(^-^ I dqe-^''^' \q)^ ® -^(|0)^ - |1),). (14) 

The oracle operates and, as seen in the previous section, 
the resulting state gains a phase (—1)-'''^*^ 

(^-^ I dq e-^^'^ (-l)/(«) ® H |1)^ . (15) 

Finally, Alice implements the final Fourier gate in her 
register (Bob's register is left unchanged and thus we 
drop it), producing the state 

|S) = ^ J dq dq' e-^'9'-2»99' (-l)/('?) \q') , (16) 

right before the measurement stage. Alice wants to mea- 
sure her initial state, so she makes a projective measure- 
ment with a little spread As around s 

ps+As/2 

fiAs= / dt\t){t\, (17) 

Js-As/2 

obtaining the state 
Ha. |S) = - r'^^^'dt I d9(-l)^(^)e-(^^+*^)'^ \t) . 

Js-As/2 J 

(18) 

Again, if the function is constant, the term (— l)-^^"?) be- 
comes ±1, leaving another Gaussian integral to solve, 
yielding 

(19) 
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after writing \t) in terms of the quadrature states \q) as 
in ( |l2| ) . The integral in t cannot be exactly solved, but as 
we are interested in states with high levels of squeezing 
(s tending to zero), the term in parenthesis can be ap- 
proximated to 1 + O(s^), the terms of order 2 and higher 
being too small to be considered. The integral in t is now 




e Ik), (20) 



which can be solved in terms of Gauss error functions 
[13 [H], finally yielding 



Has IS) = ± 



1 



\/47r 



dqT,,,{q,s)\q), (21) 



in the limit of high squeezing. The error function 
rorr('?, s) is given by 



Terr{q,s) = - <^ erf 



s - As/2 



— erf 



As/2 



(22) 

These Gauss error functions arise naturally, measuring 
the quality of the Gaussian signal in terms of the quality 
of squeezing associated to s. If the function is constant 
indeed, we should measure the initial state with total cer- 
tainty in the case of infinite squeezing as in ([6]), but due 
to the realistic Gaussian states used in practice, the prob- 
ability of measuring the initial state is given by | (s| S) p, 
and then solv ing the resulting Gaussian integrals using 



equation (16) we obtain a probability of |. This means 



that the determinism of the algorithm is lost when using 
realistic states and the algorithm becomes probabilistic, 
when the function is constant we have a 50% chance of 
getting the correct measurement at the output, we are 
no longer certain of the nature of the function with just 
one evaluation of the algorithm. However, its exponen- 
tial speedup over a classical algorithm still holds in this 
case. 



quantum algorithms. We explored the combined action 
of discrete qubits and CV states of light to perform the 
well known Deutsch-Jozsa algorithm, taking advantage 
of both schemes of computation and using only optical 
devices. In previous implementations of the algorithm, 
the main problem lies in the difficulty of producing the 
oracle operation. The hybrid model proposed here allevi- 
ates the problem by using classical feed-forward to create 
an ancillary state that performs the desired operation. 

As we have emphasized before, the importance of this 
particular algorithm lies in its relative simplicity that en- 
ables us to implement it in many settings. The hybrid 
proposal that we present also tries to inspire the imple- 
mentation of algorithms based on the Fourier transform, 
with a nontrivial oracle, in order to find the possible 
new capabilities that such algorithms would acquire in 
contrast to the available implementations. For example, 
Shor's algorithm for integer factorization satisfies both 
conditions for the optical hybrid model. 

In addition, we showed how the performance of the al- 
gorithm is diminished by considering (highly squeezed) 
Gaussian states in Alice's register. Hence, in a realis- 
tic implementation, apart from measurement errors, we 
would not have the ideal evolution of the quadrature 
states \q) either. Such implementation should provide 
some form of error correction together with repeated ap- 
plications of the algorithm for a successful, fault-tolerant 
computation that does not depend highly on the quality 
of squeezing. We sense that a better measurement stage 
would help to solve the problem. 

On the other hand, other than all-optical implemen- 
tations of hybrid computation exist, as the interaction 
of trapped atoms with a cavity mode of electromagnetic 
field [12] . These proposals may be incorporated together 
in a "super-hybrid" computer, and take advantage of all 
the schemes to achieve scalability in the computation. 



V. CONCLUSIONS 



Optical quantum computation has been proven a wide 
field of possibilities for practical implementations of 
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